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Abstract 

In a recent paper (Phys. Rev. B 78, 075316 (2008)), Sajeev and Moiseyev demonstrated that the 



j^ ' bound-to-resonant transitions and hfetimes of autoionizing states in spherical quantum dots can be 

^' 
' ■ controlled by varying the confinment strength. In the present paper, we report that such control 

(p . 

can in some cases be compromised by the presence of Coulomb impurities. It is demonstrated 

that a screened Coulomb impurity placed in the vicinity of the dot center can lead to bound- 
to-resonant transitions and to avoided crossings-like behavior when the screening of the impurity 

^ . charge is varied. It is argued that these properties also can have impact on electron transport 

O , 

Q , through quantum dot arrays. 



PACS numbers: 32.80.Zb, 73.21.La 



I. INTRODUCTION 

Autoionizing states have been very thoroughly investigated in conventional atomic sys- 
tems. The accuracy achieved nowadays in experiment and theory is as high as to resolve 
hyperfine splittings in dielectronic resonancesi. Another, currently very active, research 
topic involving autoionizing states are pump-and-probe experiments with short laser pulses, 
which make it possible to follow the autoionization process in real time and thus to resolve 
electron dynamics in atoms^''^. However, the autoionization process in artificial atoms like 
quantum dots is less known; From an experimental point of view, it is much more difficult to 
observe autoionization, since, contrary to "usual" atomic systems, quantum dots are imbed- 
ded into the semiconductor. The conventional atomic approach where autoionizing states 
are revealed as resonances in the cross section spectra for ionization by photon or particle 
impact is, in the case of quantum dots, not easily employed. Also theoretical data is rather 
scarce, since the confinment of the electrons in the dot is often modelled by a harmonic 
potential, which automatically excludes the possibility of resonances, while in reality the 
confining potential is, of course, finite. To circumvent this problem, other confinment mod- 
els were suggested to study the autoionization process in two-electron quantum dots, like a 
finite well-"^ or a Gaussian potential^i^, which indeed led to interesting observations such as 
resonance-induced enhancement of the dot sensitivity to photons^ or entanglement in reso- 
nances^. In the present work, we aim to study the effects of screened Coulomb impurities on 
the positions and lifetimes of such autoionizing resonances. The impact of charged impurities 
on the properties of quantum dots was addressed in connection with different confinment 
models, like a parabolic potential^"— or an infinite well^- but also finite potentials^"—. An- 
other interesting aspect is the behavior of quantum dots with impurities in external fields, 
which can give rise to effects like emergence of stable non- dispersive electron wave packets^, 
field-enhanced electron localization^" or change of optical properties^. However, we are 
only aware of one paper explicitly investigating the role of impurities in the autoionization 
process^. Therein, Buczko and Bassani considered a finite potential well with a hydrogenic 
impurity and chose an analytical method based on scattering theory techniques. Here, we 
adopt a different approach commonly known as "complex scaling" or "complex coordinate 
rotation" to determine the positions and widths of the resonances in a Gaussian-shaped 
spherical two-electron quantum dot with a Coulomb impurity. The method is described 



in more detail in the following section |TT] along with the computational procedure, and the 
results are presented in section lllli A discussion and conclusions are given in section IIV[ 

II. METHOD 

The Hamiltonian of a spherical Gaussian two-electron quantum dot with a Coulomb 
impurity reads 



H^-Tr^ + U„e-'-n-!^)+ f f, (1) 



.^^ ^ 2m* ^TreoeqiJ 47reoe jq^ - Qsj 

where qi,q2 are the coordinates of the two electrons, m* the effective electron mass, e 
the dielectric constant of the semiconductor, Uq the depth of the confining potential, a a 
parameter describing the range of the latter and rj the effective charge of the impurity. For 
convenience, we introduce scaled parameters as suggested in Ref.-: 

m* 
^ - -Qi, 



Vo = -^Uo, 2 

m* 

[m*) 
where rUe is the electron mass, so that the Hamiltonian can be written as: 



H=^ 



ruee^ 



7 -ir^ + Voe-^"^ + — ^ + , r 

^ V 2me 47reori/ AixeQ ri - r2 



(3) 



From here on, all quantities will be given in effective atomic units, and the values for the 
corresponding semiconductor can be reobtained from Eqs. ([2]). For example, in the case of 
GaAs we have e = 12.4 and m* /uif. = 0.067, which gives the effective energy unit 1 Ha* ~ 
11.857meV and the effective length unit Oq ^ 9.794 nm. To determine the bound as well as 
the resonant states of the system, we diagonalize the complex scaled Hamiltonian within a 
B-Spline basis set. This method was earlier applied to describe autoionizing resonances in 
atomio2^>^ and exotio^i^^ systems. The Hamiltonian (Eq.Q) is dilation-analytic, so that 
the uniform complex scaling^I"^ of the coordinates can be imposed: 

r. -^ r,e'^ (4) 



with a real parameter < ^ < 7r/4. The eigenvalues of the scaled Hamiltonian are complex. 
However, the physical states are invariant with respect to 9 while pseudo-continuum states 
are rotated into the complex plane by the angle 29. After diagonalizing the scaled Hamilto- 
nian for different values of 9, one obtains the bound states of the system as the eigenvalues 
with vanishing imaginary parts and the resonant states as complex eigenvalues, 

r 

the real part of which gives the position of the resonance and the negative imaginary part 
the halfwidth of the latter. It is connected to the lifetime r as r = h/T. The diagonal- 
ization of the scaled Hamiltonian is carried out numerically using LAPACK-routines, and 
below we briefly describe the computational procedure. The radial part of the one-particle 
Hamiltonian h is given by 

h = ^H- (-^^ + ^-(^- + 1)^' _ Voe-^^' - -^] (6) 

mge^ \ 2me drf 2merf 47reor « / 

and in the first step, the corresponding radial Schrodinger equation is solved in a sufficiently 

large box using piecewise poljTiomial functions Bk (so called B-Splines^"^ of a certain order 

k) defined on a given knot sequence. Thus, the radial one-particle eigenfunctions ipj are 

obtained in the form 

vA't) = ^CkjBk{r) (7) 

k 

with expansion coefficients Ckj- Subsequently, the matrix elements of the interaction poten- 
tial 



, V,2 = -. , 1 (8) 

m 47reo |ri — r2| 

between the obtained states are computed using multipole expansion, where the radial in- 
tegration can be carried out to machine accuracy using Gaussian quadrature, while the 
angular integration is performed analytically using Racah algebrai^^. The full Hamiltonian 
is then set up in the basis of coupled eigenstates to the one-particle Hamiltonians under 
consideration of the Pauli principle. Complex scaling is imposed by multiplying the kinetic 
terms by exp(— 2i^) and the Coulomb-terms by exp(— i^). The real and imaginary part of 
the complex rotated Gaussian confining potential are given by 



Re 
Im 



-1/oexp('-/3(rie*'')^) = -V^e^^ {-Pr'^cos{29)) cos{l3r1?>m{29)) (9) 

-1/o exp (-P [rie'^^f^ = Vo exp {~/3r^ cos(2^)) sin {/3r^ sm{29)) (10) 



and are scaled accordingly. After the setup of the Hamiltonian matrix is completed, the 
latter is diagonalized in the final step. 

III. RESULTS 

In the present calculations, we considered both donor and acceptor impurities, and there- 
fore we varied the effective charge in the regions rj > and r/ < 0. This is a quite simple 
picture of the screening mechanism, but it should be sufficient to illustrate the effect of 
impurities on autoionizing states. In a more advanced approach, one could also include a 
possible spatial dependence of the screening. For example, Kwon^i recently presented a 
model where the screening is modelled by an exponentially decreasing potential 

'^W-i^-p(-^). (11) 

where q is the true impurity charge and r^ the screening length which depends on the doping 
concentration in the semiconductor and the temperature (see Eqs.(9)-(15) in Ref.— ). We 
simplify the treatment by varying t] independently of rj, as it is often done in atomic many- 
electron systems to model the screening effects by core electrons. Also, since we are mostly 
interested in autoionizing states which are situated not very close to the dot center, this 
approximation seems reasonable. 

The parameters of the confining potential are chosen in the same region as suggested in 
Ref.-: throughout the calculations, the potential depth is kept fixed at Vq = 3 Ha* and for 
the range parameter /3 we take certain values which are well suited to illustrate the physical 
behavior we aim to demostrate. To represent the basis states, a sequence of 48 knot points 
with a box size of -R = 24 Oq is used to generate the B-Spline set, the order of which is 
k = 7 throughout the paper. We restrict our treatment to singlett resonances and include 
all configurations oi s — s,p — p and d — d type for the ^S symmetry and all configurations of 
s—p and p—d type for the ^P symmetry, which is enough to reach sufficient convergence. The 
numerical stability of the method was conffimed by successfully reproducing the positions 
and widths of the resonances as given in Ref.® for the case 77 = obtained with a Gaussian 
basis set. Figures [1] and [2] show the positions and halfwidths of the three lowest states and 
the Is-threshold for donor and acceptor impurities with different screening strengths. We 
observe that a donor impurity can turn an autoionizing state into a bound state, while for an 
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FIG. 1: (Color online) Donor impurity {rj > 0). Upper panels: energy positions of the three lowest 
states (colored lines) and the position of the Is-threshold (solid black line) as the function of the 
screening strength, shown for ^S (left) and ^P (right) symmetry. Lower panels: Halfwidths of the 
occuring resonant states. As the position of the second state (red line) crosses the threshold, it 
becomes a bound state and its width vanishes. The values for the range parameter were taken as 
/3 = 0.21 (oq)"^ for the ^S symmetry and /3 = 0.13 (oq)"^ f°^ the ^P symmetry. The energies and 
halfwidths are given both in scaled Hartree units (left axis) and meV (right axis) with material 
parameters of GaAs. 



acceptor impurity the opposite behavior is seen. Somewhat similar situations are known to 
occur in atomic systems, when a certain electronic configuration gives or does not give rise 
to autoionizing resonances depending on the nuclear charge. For example, the (ls^2p^)^5' 
state is autoionizing in beryllium^, while in the case of any other beryllium-like ion (e.g. 
beryllium-like carbon-^S) this state is bound. Of course, this analogy does not fully hold in 
the case studied here, since the confining potential remains unchanged, but to some extent 
it allows a qualitative insight in the observed behavior. Another interesting aspect in the 
case of donor impurities is that the width of the state which remains autoionizing (blue 
curve in Fig. [1]) is affected differently for different total angular momenta: in the case of 
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FIG. 2: (Color online) Same as in Fig. [T] for an acceptor impurity {rj < 0). Here, we observe 
that two states (blue dotted and red solid lines) cross the threshold. However, in the case of ^S 
symmetry, the width of the second state (red solid line) is orders of magnitude smaller than the 
one of the third state and is not shown in the plot since it is not distinguishable from zero on the 
given scale. The values for the range parameter were taken as 13 = 0.10 (flg)^^ ^^^ ^^^ ^^ symmetry 
and /3 = 0.08 (ag)"^ ^^^ ^^^ ^-^ symmetry. 

^S symmetry it is slightly increasing while for the ^P symmetry it is slightly decreasing as 
the effective impurity charge grows. Furthermore, we would like to point out that, for an 
acceptor impurity, the positions of the second and third states become very close to each 
other after crossing the threshold (blue and red curves in the upper left panel in Fig. |2]). 
By "zooming in" into the relevant region, we see that these curves show an avoided-crossing 
like behavior. It is illustrated in Fig. |3l where we also plot the energy difference AE of the 
states vs. the effective impurity charge. To summarize, we observe that the dot spectra 
are quite sensitive with respect to impurities, both concerning the positions and, in case of 
autoionizing states, also the lifetimes, in particular because impurities can cause threshold- 
crossings so that bound states become resonant or vice versa. In the following section IIVI 
we discuss the physical implications of the observed behavior which allows us to draw some 
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FIG. 3: (Color online) Left panel: "zoom-in" into the plot shown in the upper left panel of Fig. [2l 
Right panel: energy difference between the two states. The behavior is similiar to an avoided 
crossing. The range parameter of the confinment is ^ = 0.10 (ao)~^ 

conclusions. 



IV. DISCUSSION AND CONCLUSIONS 

In the context of the presented results, we focus on two main topics in our discussion: 
resonance-enhanced sensitivity of the quantum dots to photons and the role of autoioniz- 
ing resonances in transport processes through quantum dot chains. The practical purpose 
behind controlling the positions and lifetimes of resonances in quantum dots by adjusting 
the confining potential^ was a possible application of the latter as sensitive photodetectors. 
As demonstrated therein, the presence of an autoionizing resonance leads, in fact, to a very 
significant increase of the photoionization rate, since autoionizing states can be interme- 
diately populated in the photoionization process. Our results, however, indicate that in 
case of such an application attention should be paid to the purity of the semiconductor, 
since a donor impurity could turn an autoionizing state into a bound state which would 
considerably decrease the detector sensitivity. In other words, donor impurities counteract 
the controlled efficiency of such photodetectors. As for the role of acceptor impurities, we 
would like to mention their possible impact on electron transport through quantum dot 
chains. Let us for example imagine an array of quantum dots, prepared in a way that each 
of them initially contains one electron and consider the propagation of an electronic wave 
packet from one end of the chain to another—. If, in such a situation, acceptor impurities 
are present in one or several dots and give rise to autoionizing resonances, it could lead to 



8 



an additional channel for quantum transport where the electron is captured into a resonant 
state and remains there for a time span comparable to the lifetime of the resonance before 
it is released back to the continuum. This mechanism would thus compete with tunneling 
bewtween coupled quatum dots, possibly even giving rise to interference effects among the 
propagation paths. Qualitatively, one may even compare the situation to the first step in the 
process of dielectronic recombination in ions, when free electrons are captured into doubly 
excited states by simultaneous excitation of a core electron. Of course, in a semiconductor 
the resonant charecterstic would be less pronounced since the propagating electrons are not 
monochromatic; nevertheless, such a parallel between "usual" and artificial atoms is quite 
intriguing. 

In conclusion, we studied autoionizing resonances in the presence of Coulomb impurities 
in spherical Gaussian-shaped two-electron quantum dots using the complex scaled direct di- 
agonalization method. We found that donor impurities can turn resonant states with a finite 
lifetime into bound states, while acceptor impurities have the opposite effect. Implications 
of these features were discussed in the context of photoionization and transport processes in 
quantum dots, underlining the importance of the semiconductor purity in these particular 
applications. 
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